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TWO MASSIVE, LOW-LUMINOSITY CORES TOWARD INFRARED DARK CLOUDS 
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ABSTRACT 

This article presents high-resolution intcrfcrometric mosaics in the 850 /im waveband of two massive, 
quiescent infrared dark clouds. The two clouds were chosen based on their likelihood to represent 
environments preceding the formation of massive stars. The brightest compact sources detected in 
each cloud have masses w 110 M and w 60 M with radii < 0.1 pc, implying mean densities of 
(n) w 10 6 cm~ 3 and (N) w lgcm -2 . Supplementary data show these cores to be cold and inactive. 
Low upper limits to their bolometric luminosities and temperatures place them at a very early stage 
of evolution while current models of massive star formation suggest they have the potential to form 
massive stars. 

Subject headings: stars: formation — ISM: clouds — ISM: structure 



1. INTRODUCTION 

It is fairly well established that low-mass stars form 
from the gravitational collapse of dense condensations, or 
cores 1 , within molecular clouds ([Myers fc Benson 1983 



Beichman et al. 1986 



Motte et al.||1998p , and there are 
many examples ot cores in nearby star-forming regions 
(< 300 pc) that appear to represent the direct progeni- 
tors of low-mass s tars, at least in a statistical sense (e.g. 



Alves et al. 20071. However, surveys of distant regions 
where massive stars form suffer from poor physical res- 
olution, blending togeth er regions with sizes comparabl e 
to clusters of stars (e.g., Shirley et al.|2003||Evans|2008[|. 



High-resolution interferometric observations can resolve 
these massive clumps into substructure more cl osely re- 
lated to individual stars or stellar systems (e .g. , |Beuther 



et al.||2007| |Pillai et al.j|2006j jMolinari et al.||2002p . But 
in virtually all cases, these high-resolution studies have 



targeted regions with indications for young massive stars 
such as bright mid-infrared emission, masers, outflow, or 
compact Hii. 

This article presents two i nterferometric mosa i cs of in- 
frared dark clouds (IRDCs; Simon et al. 2006a Perault 

1998p selected to lack massive 



et al. 1996 Egan et al 



protostars yet likely to be pre-cluster clouds. Several 
dense cores are detected within each IRDC, but it is the 
most massive cores in each cloud are the focus of this 
article. Following a description of our observations in 
§|2] the general properties of the IRDCs are described in 
§[3] using new single dish data together with published 
data. The interferometric mosaics are then considered in 
§[4]where the physical properties of the massive cores are 
presented and discussed. 



Simon et al.| (|200"6b ) using Spitzer Galactic plan e 



vey data (|Behjamm et al. 2003 Carey et al.||2009 



and SCUBA archival data ( |jJi ErancescoeT°aT.||2do8| ) to 
be massive, dense, and lack significant 24 fira emission. 
These attributes were desired to maximize the probabil- 
ity that the clouds are in a state preceding the formation 
of massive stars. The observations are summarized in 
Table [2 

The target IRDCs were observed with the Submillime- 
ter Array 2 (SMA) in its compact configuration using the 
345 GHz receivers. Mosaic observations were designed 
to cover the brightest SCUBA emission in each IRDC. 
These observations preserved the full 4 GHz bandwidth 
and specified both high and low resolution correlator 
chunks with 0.17 and 2.8 kms -1 channel widths respec- 
tively. The data were calibrated using the MIR reduc- 
tion package 3 . Visibility phases were corrected using 
quasar observations conducted in 25 minute intervals. 
The bandpass was calibrated with observations of the 
quasar 3C279, and Uranus was used to tie down the flux 
scale to an accuracy of ~ 15%. The visibility data were 
then output into MIRIAD 4 format and inverted into the 
image domain. The mean rms level across the mosaic 
maps for MSVGOSO.SS+OO.IS and MSXG028. 53-00.25 
are 5.4 mJy and 4.6 mJy for synthesized beam sizes of 
1.9" x 1.8" and 2.0" x 1.1" oriented at position angles of 
66° and 68°, respectively. 

The James Clerk Maxwell Telescope 5 (JCMT) 
equipped with the HARP-B receiver array and ACSIS 
backend ( Dent et al.||2000[ ) was used to observe our 
target IRDCs over several nights. Spectral lines of CO, 
13 CO, C 18 0, N 2 H+, and H 2 D+ were placed in spectral 
windows of 250 MHz width and 61 kHz resolution. Flux 



2. OBSERVATIONS AND DATA REDUCTION 

The infrared dark clouds MSXG030.88+00.13 and 
MSXG028. 53-00.25 were chosen from the catalog of 

Electronic address: js@ifa.hawaii.edu 

1 Following the convention presented in jWilliams e t al.| ( |2000 1 , 
"cores" describe dense condensations of molecular gas with M < 
100 Mq that will form individual stars or small stellar systems, 
while "clumps" describe dense cloud structures with mass > 
1000 Mq that are more likely to form clusters of stars 



2 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 

3 |http : //www, cf a .harvard ■ edu/~cqi/mircook.html| 

4 http : //bima. astro .umd. edu/miriad 

5 The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facil- 
ities Council of the United Kingdom, the Netherlands Organisa- 
tion for Scientific Research, and the National Research Council of 
Canada. 
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TABLE 1 

Summary of Observations 



Date Waveband Conditions Target a Calibrators Obs. Type 



The Submillimctcr Array 



2007Jun24 


350.9 GHz b 


T225 


SJ 0.08 c 


B 


1751+096/1743-038 


1 


2007Jul01 


351.0 GHz b 


T225 


SS 0.08 c 


B 


1751+096/1743-038 


1 


2007Jul07 


350.9 GHz b 


T225 


« 0.08 c 


A 


1751+096/1743-038 


1 


2007Octl9 


350.9 GHz b 


T"225 


ss 0.08 c 


A 


1751+096/1743-038 


1 


2008Jun02 


340.5 GHz b 


T225 


» 0.15 c 


A,B 


1751+096/1911-201 


1.2 


2008Junl5 


340.5 GHz b 


T"225 


» 0.11 c 


A,B 


1751+096/1911-201 


1,2 



James Clerk Maxwell Telescope 



2007Jull7 


334.8 GHz/350.7 GHz b 


T225 


« 0.10 c 


A,B 


V437Sct/16293-2422/G45.1 


3 


2007Jull8 


350.7 GHz/367.3 GHz b 


T"225 


« 0.06 c 


B 


V437Sct/16293-2422/HD235858 


3,4 


2007Jull9 


334.8 GHz/350.7 GHz/367.3 GHz b 


T"225 


« 0.06 c 


A,B 


V437Sct/16293-2422/G34.2 


3,4 


2007Jul20 


334.8 GHz/350.7 GHz/367.3 GHz b 


T225 


ai 0.07 c 


A,B 


V437Sct/16293-2422 /G34.2 


3,4 


2008Aug06 


367.3 GHz b 


T225 


« 0.05 c 


A 


WAql 


4 


2008Novl2 


367.3 GHz b 


T225 


sj 0.07 c 


A 


HD179821/NMLCyg 


4 



Canda France Hawaii Telescope 

2008AuglO J, H, K s S « 0.6" d A^B 2MASS 5 

2008Augl6 J, iC s 6 a sj 0.7" d A,B 2MASS 5 

Note. — (1) Mosaic observations, (2) Targeted (single pointing) observations, (3) raster scan mapping, (4) jiggle chop 
observations (5) direct imaging. 

a Target A corresponds to M5AG030. 88+00. 13, and target B corresponds to M5AG028. 53-00.25. 

b Local oscillator frequency. 

c Atmospheric opacity at 225 GHz. 

d Full width at half maximum of the point spread function. 



levels are expected to be accurate to 15% and the 
pointing accuracy of the final maps are estimated to 
be better than 3.5". The Starlink software suite was 
used to grid and output the data into FITS format 
and IDL was used to perform final calibrations and 
coadding. A main beam efficiency 77MB = 0.7 is used to 
convert corrected antenna temperatures to main beam 
temperatures. 

The Canada-France Hawaii Telescope 6 (CFHT) was 
used to observe o ur target clouds u sing the WIRCam 
infrared detector (Puget et al. 20041. The observations 



were carried out in queue mode over the course of two 
photometric nights. The pre-processed data (de-biased, 
flat-fielded and sky-subtracted by the CFHT pipeline) 
were downloaded and further reduced using the TER- 
APIX 7 software suite. The images were registered, com- 
bined by weighted mean, and tied to the 2MASS 8 point 
source catalog flux scale to an accuracy of better than 
0.02 magnitudes. 

3. THE IRDC ENVIRONMENTS 

Figures |TJl and [TJ) show images from the Spitzer 
GLIMPSE survey in which the mid-infrared extinction 
features that define the clouds can be clearly seen. Con- 
tours of velocity integrated 13 CO (3 — 2) overlay the im- 
ages in yellow. M5XG030. 88+00. 13 has 2 velocity com- 

6 The Canada-France-Hawaii Telescope (CFHT) is operated by 
the National Research Council of Canada, the Institut National 
des Sciences de l'Univers of the Centre National de la Recherche 
Scientifique of France, and the University of Hawaii. 



http : //terapix . iap . f r 

1'he Two Micron All Sky Survey is a joint project of the Uni- 
versity of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 



ponents along the line-of-sight with centroids at 95 and 
107 km s -1 with respect to the local standard of rest. 
The higher velocity component shown in Figure TJi at 
13 CO emission levels of (5 + 3 n) Kkms -1 (n — 0, 1,2...) 
is widespread across the region suggesting that the IRDC 
is part of a larger complex including the bright infrared 
cluster to the northwest. MSXG028.53-00.25 has a 
single velocity component in 13 CO at 87kms~ 1 shown 
at contour levels of (6 + 2n)Kkms _1 seen to be spa- 
tially confined to the region of mid-infrared extinction. 
The presence of high volume density gas in these clouds 
is confirmed with detections of N2H + (4 — 3) shown as 
magenta contours at levels of (1 + 0.5n)Kkms _1 and 
(1 + 0.4 n) Kkms" 1 for Figures [la andflp, respectively. 

The velocity centroid of M5XG030. 88+00. 13 places it 
at the tangent point of Galactic rotation in this direction. 
Therefore a distance of 7.2 kpc can be derived using the 
simple geometric relationshi p d = Rg cos (7) cos (6), where 



Ro is taken to be 8.4 kpc (Ghez et al. 
20091. The distance to MS 



Reid et al 
is taken to 

be 5.4 kpc (iRathborne et al. 2006 1 . The errors on these 



kinematic distances may be 15% o r more due to non 
circul ar motions in the Galaxy (e.g., Roman-Duval et al. 
2009 1 . The masses of the IRDCs are estimated to be 



several thousand solar masses based on CO isotopologue 
emission assuming 15 K gas in local thermodynamic equi- 
librium, SCUBA 850 /im continuum emission, and virial 
equilibrium. 

There is no A20 cm or A6 cm emission detected tow ard 
either IRDC rtHelfand et al.||2006| [White et al.| p005], A 



total A20 cm flux ot 250 mJy seen toward the 1R bright 
cluster in Figure [lk suggests the presence of a B0 star 
(Giveon et al. 2005b. Two class II methanol masers are 
seen in the field of MSXG030. 88+00. 13 with velocities in 



rough agreement with the cloud ( jPestalozzi et al. 2005 ) 



Two Massive, Low-Luminosity Cores 
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MSX GQ30. 88 + 00.13 M SX G028.53-00.25 




18 h 47"'18 s 15 s 12' 18 h 47 m 09" 18 h 44 m 21* 18' 15* 18" 44™ 12" 

Right Ascension (J2000) Right Ascension (J2000) 



Fig. 1. — Gray scale images of 5.8 + 8.0 (im emission toward two infrared dark clouds. Contours of 13 CO(3 — 2) emission are shown in 
yellow and contours of N2H + (4 — 3) are shown in magenta. The circled crosses d enot e the positions of methanol (white) and water (cyan) 
masers. The positions of the brightest compact cores detected by the SMA (see §|4.1|l are shown with the green cross. 



No class II masers are seen in the MSXG028. 53-00. 25 
field. However, one water maser has been detected to- 
ward an embedded source in this cloud ( |Wang et al. 
2006 1 . Fits to the spectral energy distributions of MIPS 
24 /mi sources within the IRDCs suggest t hey are not 



massive stars (Robitaille et al. 2006 20071. Although 



given the large column depths observed toward the 
IRDCs and their distances, it is possible that more pro- 
tostars are embedded in the cloud than are detected in 
the Spitzer data. We estimate that less than ~ 100 stars 
are embedded within the IRDC based on extrapolations 
from mid-infrared emission from nearby, intermediate 



mass star forming regions (Padgett et al. 2008 Rebull 



et al. 20071 and a S alpeter stellar mass funct ion w ith a 



lognormai turnover ( |Salpeter||1955| |Chabrier||2003[ ) 

4. HIGH-RESOLUTION SUB-MILLIMETER WAVE 
IMAGING 

Figures J2a and[2]D show the contours of the SMA mo- 
saics overlaid on Bpitzer mid-infrared composite images. 
The brightest sources in each mosaic dominate the com- 
pact sub- millimeter wave emission toward the IRDCs. 
These massive cores lie near the near the peak of the 
SCUBA emission and are isolated with respect to other 
compact sources. They are also isolated in mass with the 
next most massive cores in each cloud being more than 
a factor of 3 less massive. Table [2] displays the physical 
characteristics of the two compact submillimeter sources. 



4.1. Two Massive, Low-Luminosity Cores 

The derived masses of the two cores are « 110 M@ 
and w 60 M Q with densities of (n) « 10 6 cm -3 and 
(N) « lgcm~ 2 . Their free fall, or dynamical timescales 
are thus a few 10 4 yrs implying that they are most likely 
short lived structures. Both cores appear marginally re- 
solved, but are detected in the longest baseline data sug- 
gesting that perhaps there are resolved and unresolved 
components. 

Figure [3] shows the spectral energy distributions of 
the cores consisting of upper limits from 1.25 /mi to 
70 /mi and 853 /mi fluxes from Table [2| A best fit gray 



body curve with r = (v/v c ) 2 , and v c = 6THz ( |Ward- 
|Thompson et al.||2002 1 limits the envelope temperatures 
while summation under the solid curves of Figure [3] pro- 
vide upper limits to the bolomet ric luminosities ana tem- 
peratures ( Myers fc Laddfl993 l 



High values of M em /L^ x « 2.8 in solar units and 
low values for Tboi < 30K for these cores are signs of 
extreme youth. Class protostars, t he earliest SEP 
class, are marked by M e nv /L®-® > 0.4 ( Bontemps et al. 
19961 lAndre et al.||2000| and T hol < 70 K (jChen et al. 



1995). Given the large distances to these cores and the 



photometric sensitivity, the existence of a low luminos- 
ity sources deeply embedded within the cores cannot be 
ruled out. However, it is clear that no massive protostars 
exist in these cores. 

Across the SMA bandpass, the cores are only detected 
in CO(3-2) and HCO+(4-3). The CO(3-2) is strongly 
filtered by the interferometer making a useful interpre- 
tation of the emission from this low density gas tracer 
difficult. However, no clear indication of outflow is seen 
toward either core in the single dish or interferometric 
data further indicating their extreme youth. Upper lim- 
its to SO(8 9 - 7 8 ) and CH 3 OH(13 7i 7 - 12 7 , 6 ) emission 
highlight the difference between these cores and t he spec- 
tral s ignatu res of hot cores such as Orion KL ( |Schilke| 



et al. 1997 



a- 



Figure 14] shows the composite spectra of molecular 
transitions tracing dense gas toward the two cores with 
the velocity scale shifted relative to the centroid velocity 
of N2H + (4— 3) emission measured to be 107.2kms _1 and 
86.8 km s -1 for core 1 and 2, respectively. Core 1 shows 
an HCO + emission feature that lies blueward of the sys- 
temic velocity. This may be an indication of HCO + self- 
absorption, but there is also the possibility that the emis- 
sion in the central channels is widespread and filtered out 
by the interferometer. Core 2 shows a weak feature in 
the HCO + spectral window blueward of its systemic ve- 
locity that may be a sign of a similar spectral feature to 
core 1. 

Also shown in Figure [4] are spectra of ortho- 
H2D + (li.o — li.i) toward the cores. Core 1 shows 
an emission feature with 4.3 a significance centered at 
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MSX G030.8B + 00.13 Msx G028.53-00.25 




18 h 47 m 18 !l 15* 12 s 18 h 47 m 09 s lB 1 '44 m 21 s ' 16 s 15 s lB h 44 m 12 s 

Right Ascension (J2000) Ri S ht Ascension (J2000) 



Fig. 2. — Spitzer composite images (red = MIPS 24 (im, green = IRAC 8.0 fim + 5.8 ^m, blue = IRAC 4.5 fim + 3.6 fim) overlaid with 
contours of sub- millimeter continuum emission from the SMA mosaic data at levels of 3a n (orange) and SCUBA archival data shown at 
levels of 250 n mjy (gray) with n = 1, 2, 3, ... The coverage of the SMA and SCUBA data are outlined with a dashed and dotted contour, 
respectively. 



TABLE 2 
Massive, Low-Luminosity Cores 





R.A. 


Decl. 






Mass a 


r b 


<n> 


(N) 




Tbol 


Lhoi 


Name 


(J2000) 


(J2000) 


(mjy/bm) 


(mjy) 


(M ) 


(pc) 


(10 6 cm" 3 ) 


(gem 2 ) 


(K) 


(K) 




Core 1 


18:47:13.7 


-01:45:03.7 


119 


255 


110 


0.082 


0.82 


1.1 


< 19 


< 31 


< 460 


Core 2 


18:44:18.0 


-03:59:23.0 


118 


224 


(50 


0.056 


1.5 


1.3 


< 18 


< 30 


< 170 



|1994| Table 1). 

b Radii are derived by deconvolving the synthesized beam from the effective radius, r = >J R?g — fl^ eam , where R c g 
and A is the area contained within a 3 it contour. 



lOe.ekms- 1 with a full width of 0.9 ± 0.3 kms" 1 . Core 
2 shows a broader feature with 2.0 ± 0.7kms _1 width 
at a 3.1 a significance level. This molecular species is 



in the deei 


3 interior of pre- or proto-stellar cores 1 


Walm- 


sley et al. 


2004 


Vastel et al. 


2006 


Caselli et al. 


2008}. 



T bol <31K 
r bo , < 30 K 




A (/im) 



Fig. 3. — Spectral energy distributions of core 1 (black) and core 
2 (red). The detections at 853 fim have error bars smaller than the 
symbol size, while upper limits at wavelengths from 1.25 to 70 fim 
arc shown. 



Following |Caselli et al. (20081, the integrated line fluxes 
of 0.13Kkms _i and 0.17Kkms _1 for core 1 and 2 re- 
spectively translate to column densities of A^H^D^ 
10 1 



111 cm 2 using a kinetic temperature Tk = 15 K and a 
critical density n CT — 10 6 cm~ 3 . The fractional abun- 
dances of ortho-H 2 D+ are iV(H 2 D+)/iV(H2) « 3-5 x 
10" 13 . Varying T k from 10 to 15 K and N rr from 10 5 



are iV(H 2 D H 
Varying Tk from 10 to 15 K and N C1 
to 10 6 cm -3 changes these values by a factor of about 2. 
These numbers are significantly smaller than the values 
fo und for low-mass pre- and pr oto-stellar cores ( |Vastel| 

f aLl|2006| IC asell i et al.||2008| ) as well as in Orion ki 
arju et al.||2006p . This could be due to an intrinsic 
dearth of rl2lJ + . However, given the poor physical reso- 
lution of our single dish observations and the propensity 
for H2D + to trace the innermost regions of cores, it may 
be that this discrepancy is due in large part to beam 
dilution. 

4.2. Discussion 

The signposts of massive star formation toward 
MS'VG030.88+00.13 lend credence to the idea that there 
might exist regions within the cloud where massive stars 
will, but have not yet, formed. The characteristics of 
core 1 make it a good candidate for one such region. 
While M5XG028.53-00.25 currently shows no signs of 
high-mass star formation, core 2 shares similar traits to 
core 1 as being a potential massive star precursor. 

The bolometric luminosity of cores versus the enve- 



Two Massive, Low-Luminosity Cores 
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T3 

(L) 

O 

to 

A 
V 




-20 -10 

Velocity (km/s) 



10 



Fig. 4. — Average spectra of high-dipole moment molecules to- 
ward core 1 (top) and core 2 (bottom) in molecular transitions 
tracing high volume density gas. The N2H+ and H2D+ spectra 
are scaled by a factor of 10 for the plot. Velocities are with respect 
to the centroid of N2H + emission for each core measured from a 
Gaussian fit to the line profiles. 

lope mass, or the L^ i-M env diagram is a useful param- 
eter space to visualize the evolution of cores into stars 
(e.g., Bontemps et al. 1996). Figure [5] shows data from 
several previous studies of massive star formation plot- 
ted on the Lb i-M env diagram (see figure caption). The 
linear trend in log-log space is typical for such studies. 
Evolutionary tracks b ased on the turbulent core model of 
McKee fc T an] ( |2003| ) calibrated by the data of|Molinari' 
are also overlaid on the plot 



et al.| ( |2008p are also overlaid on the plot. The positions 
of core 1 and 2 from this study are seen to lie well below 
the trend set by previous data, another indication that 
the y are at an ear l y ev olutionary stage. According to 
the Molinari et al. ( |2008[ ) evolutionary tracks, core 1 is 
on track to evolve into a star with iboi ~ 10 4 L Q , or of 
an early B spectral type while core 2 is set to evolve into 
a star with Lboi ~ 2 x 10 3 L©, or a mid-B spectral class. 

The final stellar masses implied from the evolutionary 
tracks in Figure [5] may be underestimates. The star for- 
mation efficiency of dense cores is expected to be be- 



2ncy 
% [I 



tween 25-50% (iMatzner & McKee 2000 Alves et al. 



20071. Therefore if these cores collapse monolithically 
they have sufficient mass to form stars with mass be- 
tween 15 and 60 Mq. Indeed, the evolutionary models of 



Andre et al.| ( 2008 1 put core 1 on track to evolve into a 
star of « 50 Mq and core 2 into a star with w 25 Mq . 
However, it is not possible to tell from our data whether 
these cores will undergo monolithic collapse or will fur- 
ther fragment into smaller cores. Observations resolving 
the expected fragmentation length scales (~ 0.03 pc) are 
needed. 

The large millimeter continuum survey of the Cygnus 
X region has also reveale d dense and massi ve cores at an 
early stage of evolution ( Motte et al. |2007 1 . However, of 
the most massive cores in their sample (M > 40 Mq ) , 
61% have 21 /im MSX sources associated with them, 
and the remaining cither contain compact Hn, show 
70 fim emission above our scaled detection limits (S. Bon- 
temps, F. Motte, private communication), or have SiO 
line strengths and widths that indicate protostellar ac- 
tivity. 

Therefore it seems that the cores of this study are 
unique, and their existence may have implications for our 
understanding of massive star formation. The selection 
criteria and use of high-resolution in the sub-millimeter 
wave band are key components to these disc overies. As 
the charact erizations of IRDCs progress (e.g., Chambers 
et al.| 2009 ) further high-resolution studies of pre-cluster 
environments will provide statistics on the frequency and 
nature of these kinds of objects. 



5. SUMMARY 

This article presents intcrferometric mosaic observa- 
tions in the 850 /im waveband toward two infrared dark 
clouds chosen to be likely representations of pre-cluster 
environments in the Galaxy. Both clouds are massive and 
dense but show no sign of massive star formation. The 
most massive cores in each mosaic dominate the compact 
850 ^m emission and are spatially isolated near the peak 
of the low-resolution continuum emission. With masses 
of 110 and 60 Mq average densities of (n) ss 10 6 cm -3 
and (N) w lgcm~ 2 they have the potential to form 
massive stars. However, upper limits to their bolomet- 
ric temperatures and luminosities, no clear indication of 
CO outflow, and their relatively featureless spectra all 
indicate a very early stage of evolution. Detections of 
ortho-H 2 D+(li i0 - 1 M ) and N 2 H+(4 - 3) support the 
interpretation that these cores are in a cold and dense 
state preceding significant star formation activity. The 
comparison of these data with theoretical and observa- 
tional studies of massive star forming regions suggest 
that these cores occupy a unique region of Lboi _ -^env 
parameter space placing them on track to evolve into 
stars anywhere from mid to early B stars up to O stars. 



The author is grateful for the referee's suggestions as 
well as the helpful input provided by many people in- 
cluding David Jewitt, John Johnson, Sylvain Bontemps, 
Frederique Motte, Jonathan Williams, Emeric Le Floc'h, 
Qizhou Zhang, Steve Longmore, Thushara Pillai and Jill 
Rathborne. 
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